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ABSTRACT 

Over half of GeV gamma-ray sources observed by the EGRET experiment have not yet been identified 
as known astronomical objects. There is an isotropic component of such unidentified sources, whose 
number is about 60 in the whole sky. Here we calculate the expected number of dynamically forming 
clusters of galaxies emitting gamma-rays by high energy electrons accelerated in the shock wave when 
they form, in the framework of the standard theory of structure formation. We find that a few tens 
of such forming clusters should be detectable by EGRET and hence a considerable fraction of the 
isotropic unidentified sources can be accounted for, if about 5% of the shock energy is going into electron 
acceleration. We argue that these clusters are very difficult to detect in x-ray or optical surveys compared 
with the conventional clusters, because of their extended angular size of about ~ 1°. Hence they define 
a new population of "gamma-ray clusters" . If this hypothesis is true, the next generation gamma-ray 
telescopes such as GLAST will detect more than a few thousands of gamma-ray clusters. It would provide 
a new tracer of dynamically evolving structures in the universe, in contrast to the x-ray clusters as a 
tracer of hydrodynamically stabilized systems. We also derive the strength of magnetic field required for 
the extragalactic gamma-ray background by structure formation to extend up to 100 GeV as observed, 
that is about 10~ 5 of the shock-heated baryon energy density. 

Subject headings: cosmology: theory — diffuse radiation — large-scale structure of the universe — galaxies: 
clusters: general — gamma rays: theory 



1. INTRODUCTION 

The deepest image of the universe in the high energy 
gamma-ray band beyond 0.1 GeV has been obtained by 
the EGRET experiment (Hartman et al. 1999). Identi- 
fied sources of the third EGRET catalog include the Large 
Magellanic Cloud, five pulsars, and 66 active galactic nu- 
clei (AGNs) of the blazar class. However, over half of the 
EGRET sources have not yet been identified as known 
astronomical objects, and their origin is one of the most 
interesting mysteries in astrophysics. The distribution of 
these unidentified sources can be interpreted as the sum of 
the Galactic component along the Galactic disc (|6| < 40°) 
and another isotropic (i.e., likely extragalactic) component 
(Ozel & Thomson 1996; see also Fig. 2 of Mukherjee 
et al. 1995). Several candidates have been proposed as 
the origin of the Galactic unidentified sources, including 
molecular clouds, supernova remnants, massive stars, and 
radio-quiet pulsars (see e.g., Gehrels & Michelson 1999 
and references therein). However, almost no candidate 
has been proposed to explain the extragalactic unidentified 
sources, except for undetected AGNs. Recently, Mirabal 
et al. (2000) have performed comprehensive follow-up ob- 
servations for one of the high-latitude unidentified EGRET 
sources (3EG J1835+5918, b = 25°) in X-ray, optical, and 
radio wavebands. They found that any known class of GeV 
gamma-ray sources including blazars and pulsars cannot 



be the origin of 3EG J1835+5918, and it suggests that this 
source belongs to a new class of GeV gamma-ray emitters. 

There are 19 unidentified sources with high galactic lat- 
itude of \b\ > 45° (about 60 in the whole sky), and 7 of 
them are noted as 'em' in the third EGRET catalog, i.e., 
possibly extended or multiple sources that are inconsistent 
with single point sources. Although this 'em' designation 
is quite subjective and we should be careful in interpret- 
ing this result (see Hartman et al. 1999 for detail), this 
may indicate that there is an extended and extragalactic 
population in the unidentified EGRET sources. Recently 
Gehrels et al. (2000) chose 'steady' unidentified sources 
from the third EGRET catalog, to eliminate any source 
with high variability typical of flaring AGNs. These steady 
sources are mostly distributed at low and mid galactic lat- 
itude < 40°) and hence should be the Galactic origin, 
but 7 sources are still located at \b\ > 45°, that may be 
steady extragalactic sources whose number is ~ 24 in the 
whole sky. Therefore, steady astronomical objects with an 
extended nature are worth being investigated as a possible 
origin of unidentified EGRET sources. 

It is widely believed that the observed structures in the 
universe have been produced via gravitational instability. 
Currently the most successful theory of structure forma- 
tion is the cold dark matter (CDM) scenario, in which 
the structures grow up hierarchically from small objects 
into larger ones. When an object collapses gravitationally 



2 



and virializes, the baryonic matter in the object is heated 
by shock waves up to the virial temperature. Particles are 
expected to be accelerated to high energy by shock acceler- 
ation, and accelerated electrons scatter the photons of the 
cosmic microwave background radiation (CMB) to high 
energy gamma-ray bands by the inverse-Compton mecha- 
nism. Existence of such nonthermal electrons is inferred 
from radio and hard x-ray observations for some clusters 
of galaxies (e.g., Fusco-Femiano ct al. 1999), although the 
origin of the nonthermal electrons is not yet clear. It has 
also recently been argued that this radiation process in the 
intcrgalactic medium may explain the diffuse extragalac- 
tic gamma-ray background radiation (EGRB) observed in 
the EGRET range (Loeb & Waxman 2000). However, it 
is still highly speculative and difficult to test whether this 
process is really the origin of the EGRB, since the contri- 
bution by unresolved active galactic nuclei is also of the 
same order of magnitudes (see, e.g., Miicke & Pohl 2000 
and references therein). 

On the other hand, if the structure formation is actu- 
ally an efficient radiation process of gamma-rays, clusters 
of galaxies should be strong emitters of gamma-rays when 
they dynamically form, and the detectability of such form- 
ing clusters as discrete sources is of great interest as a new 
probe of structure formation in the universe as well as a 
test for the scenario proposed by Loeb & Waxman (2000) 
for EGRB. In this paper we make a theoretical estimate of 
the number and angular size of such gamma-ray emitting 
clusters detectable by EGRET, based on the standard the- 
ory of structure formation in the CDM universe. We find a 
few tens of such forming clusters should have already been 
detected by EGRET. Detectability of such forming clus- 
ters in other wavebands such as optical or x-ray bands will 
be discussed, in comparison with the conventional clusters 
of galaxies identified in these wavebands. We will also cal- 
culate the EGRB spectrum from structure formation and 
derive a quantitative relation between the higher cut-off 
photon energy and magnetic field strength. 

Throughout this paper, we assume a CDM universe with 
the density parameter Q — 0.3, the cosmological constant 
fi A = 0.7, the Hubble constant h = iJ /(100km/s/Mpc) = 
0.7, the baryon density parameter fig — 0.015ft. -2 , and 
the density fluctuation amplitude erg = 1- These cosmo- 
logical parameters are consistent with various observations 
including those of the CMB fluctuations (e.g. de Bernardis 
et al. 2000) and the abundance of x-ray clusters of galaxies 
(e.g. Eke, Cole & Frenk 1996; Kitayama & Suto 1996b). 

2. GAMMA-RAY LUMINOSITY AND FLUX FROM FORMING 

CLUSTERS 

We first estimate the gamma-ray flux of a gravitation- 
ally bound object of total mass M that virializes at red- 
shift z. The typical radius r v i r , density p v ; r , circular ve- 
locity V c , and temperature T v j r of the object can be com- 
puted from the spherical collapse model (Peebles 1980; 
Kitayama & Suto 1996b), that is widely used in study 
of structure formation. The total gravitational energy 



given to the baryon gas in the forming cluster is given 
by -Ebaryon ~ (3/4) (Qb/^o)MV^ . It is reasonable to ex- 
pect that a fraction £ e ~ 0.05 of this energy goes into 
accelerated electrons, since such a fraction is inferred for 
acceleration of cosmic ray electrons in a supernova rem- 
nant SN 1006 from x-ray and TeV observations (Koyama 
et al. 1995; Tanimori et al. 1998) and consistent with 
the energetics among cosmic-rays, turbulent motions, and 
supernova rate in our Galaxy. It has also been suggested 
that the diffuse radio and hard x-ray emissions observed in 
the Coma cluster (and possibly other several clusters) can 
be attributed to nonthermal electrons with the electron 
energy fraction of the same order (e.g., Fusco-Femiano et 
al. 1999). Therefore, we use £ e = 0.05 to determine the 
normalization of electron energy spectrum throughout this 
paper. 

The maximum Lorentz factor of electrons is constrained 
by the competition of the Fermi acceleration time and cool- 
ing time by inverse-Compton (IC) scattering of CMB pho- 
tons. The acceleration time is given by 

t m 

^acc y 2 \±) 
v s 

= 1.6 x IQ-^B^Vj yr, (2) 

where tl = m e ~f e / (eB) is the Larmor radius of elec- 
trons, 7 e the electron Lorentz factor, e the electron 
charge, B^q = B/(l(iG) the magnetic field, and V s ^ = 
T4/(10 3 km/s) the shock velocity that is of the same order 
of magnitudes with the circular velocity of a halo, V c . On 
the other hand, the IC cooling time is 

7 e m e c 2 

tW = (4/3)ca T ^CMB7e 2 

= 2.3 x lO 12 7" 1 (l + 0)- 4 yr, (4) 

where <7t is the Thomson cross section and E/qmb = 4.32 x 
10~ 13 (l + z) 4 erg cm~ 3 is the CMB energy density. Equat- 
ing these expressions of i acc and tic, we have the maxi- 
mum value of 7 e as 7 e , max = 1-2 x 10 8 (1 + z)~ 2 B 1 Jq V Si z. 
We assume the energy distribution of accelerated electrons 
as a power-law with an exponential cut-off at 7e jma x, i-e., 
dNe/dje oc j~ a cxp(— 7e/7e,max), with the standard par- 
ticle acceleration index of a ~ 2. As mentioned above, the 
normalization of this spectrum is set by the equation 

J dj e m e C 2 ~/ e ^- =£ e ^baryon , (5) 

with the parameter £ e = 0.05. The observed photon en- 
ergy e 7 scattered by electrons is related to 7 e as e 7 = 
(4/3)7 2 ecMB,o, where ccmb,o = 6.4 x 10~ 4 cV is the mean 
photon energy of the CMB at z — 0. 

The cooling time of electrons corresponding to pho- 
ton energy e 7 can be written as tic = 2.1 x 
10 6 (e 7 /GeV)~ 1/,2 (l+z)~ 4 yr, and this should be compared 
with the time for the shock wave to propagate the radius of 
the virialized halo, t s hock ~ r v i r /V s = r v j r /(4V c /3). Here 
we have estimated the shock velocity as V s = (4/3)14, 
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that is a velocity of a strong shock when a material is 
shocked by a supersonic piston with a velocity of V c , i.e., 
a typical bulk velocity of material in a collapsed halo. By 
using the spherical collapse model mentioned above to cal- 
culate r v ; r and V c for the halo, this time scale can be writ- 
ten as t shock = (3/4) 3 / 2 ^- 1 /2( Gpvir )-i/2 „ 1.5(1 +z )-3/2 

Gyr, that is essentially the dynamical time of the halo. 
Note that it depends only on the redshift and not on the 
halo mass. From this argument the cooling time of elec- 
trons emitting gamma-rays above 0.1 GeV is always much 
shorter than the time scale i s hock during which the shock 
is alive and a halo is an active gamma-ray emitter. Hence 
the total number of gamma-rays emitted from a forming 
halo during the time i s hock is given as: 

dN J (e J ; M, z) _ m e j e dN e dj e 
de 7 e 7 G?7 e c?e 7 ' 

and the observed photon flux of gamma-rays during the 
shock propagation time is 

dF(e 7 ; M,z) _ (1 + z) dN 7 1 
de 7 Andl de 1 t shock 

where g?l is the standard luminosity distance. 

We introduce a parameter £s to determine the mag- 
netic field, which is the fraction of magnetic energy den- 
sity in the total baryon energy density of the halo, i.e., 
B 2 /(8n) = 6?£baryo„/(47rr3 ir /3). Then by using the 
spherical collapse model again, the magnetic field of a clus- 
ter can be written as 

B ~ ai7 (^) 1,2 (jo^k)" 3(1+s)2 " G ^ ,8) 

A magnetic field of ~ 0.1-1 fiG is often observed in intr- 
acluster medium of rich clusters (Kronberg 1994; Fusco- 
Femiano et al. 1999; Rephaeli, Gruber, & Blanco 1999), 
and hence we use £b = 10~ 3 to be consistent with the ob- 
servations. It should be noted, however, that this param- 
eter is important only for the maximum photon energy of 
the gamma-ray spectrum (well beyond 10-100 GeV with 
£b = 10~ 3 , see §Q), and the gamma-ray flux above 100 
MeV measured by the EGRET is almost insensitive to 
this uncertain parameter, unless the particle index a sig- 
nificantly deviates from the standard value of 2. 

It is interesting to apply the above model to nearby 
known clusters of galaxies. For example, the Coma clus- 
ter has the total mass of M ~ 1O 15 M0 and is located 
at z = 0.023. The flux with these parameters becomes 
F(> O.lGeV) ~ 6.5 x 10" 7 photons cnT^sec -1 , that is 
about 10 times brighter than the observational upper limit 
on this cluster by EGRET, 4 x 10~ 8 photons cm _2 sec _1 
(Sreekumar et al. 1996). This does not mean, however, 
that our model is incorrect. Our model is relevant for just- 
forming clusters of galaxies in which the violent shock gen- 
erated by gravitational collapse is still alive. The Coma 
cluster is thought to have formed more than a few dynam- 
ical times ago and is now hydrodynamically stable after 



violent shock has disappeared. Gamma-ray flux can then 
be much weaker than our estimation. On the other hand, 
it suggests that there was an epoch during which this clus- 
ter was a strong gamma-ray emitter, and also that there 
may be other clusters visible by EGRET, that are just dy- 
namically forming and have not yet reached hydrostatic 
equilibrium. 

3. EXPECTED NUMBER OF GAMMA-RAY CLUSTERS 
DETECTABLE BY EGRET 

The number of such forming clusters of galaxies with 
flux stronger than F can be calculated as 

//>oo A\f 
dz dM — i? form (M, z) Uock ,(9) 

JM(z; F) dz 

where dV/dz is the comoving volume element of the uni- 
verse, i?form the formation rate of dark haloes (or clus- 
ters) per unit mass, cosmic time, and comoving volume, 
and M(z; F) the mass of a cluster collapsing at redshift z 
whose flux is F. Here we have taken into account that clus- 
ters are active gamma-ray emitters only during the time 

^shock- 

There are several formulae to calculate the formation 
rate i?f r m in the framework of the standard theory of 
structure formation. As is well known, the Press & 
Schechter (1974, PS; Peebles 1980) formalism provides a 
formula of mass function (i.e., number density of haloes as 
function of mass and redshift), that is in reasonable agree- 
ment with TV-body simulations. Here we want the forma- 
tion rate of haloes rather than the mass function at a given 
epoch, because we need to calculate the number of collaps- 
ing objects experiencing shock at each epoch. A naive pre- 
scription to obtain this quantity is to take a time derivative 
of the PS mass function, RdPS, although this is not exactly 
-Rform but rather interpreted as RdPS = ^form — -Rdest, 
where i?dcst is the rate of destruction of haloes by merging 
into even larger structures. Consequently, i?dPS becomes 
negative at small mass scales where i?dost is significant. 
As shown below, however, the number of objects visible 
by EGRET is dominated by massive clusters forming in 
recent past (z <C 1) that are the largest structures in the 
universe, and hence i?dcst is negligible. Therefore it is a 
reasonable approximation to use RdPS i n the mass range 
where it is positive. Alternatively, there are several exten- 
sions to the PS theory for computing i?f or m (e.g., Lacey & 
Cole 1993; Sasaki 1994; Kitayama & Suto 1996a; Percival 
& Miller 1999). In what follows we use RdPS and the for- 
mulae of Sasaki (1994) and Kitayama & Suto (1996a) to 
take account of theoretical uncertainties in i?f or m- 

Figure [l] shows the theoretically predicted log iV-log F 
of forming clusters. At least a few tens of clusters should 
be visible by the EGRET, and a significant fraction of the 
isotropic unidentified EGRET sources can be accounted 
for. It is also interesting to note that the predicted num- 
ber is similar to that of 'em' isotropic unidentified EGRET 
sources, i.e., possibly extended sources. The predicted 
number at the EGRET flux limit also agrees with the num- 
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10 -11 1Q -10 1Q -9 10-8 1Q -7 1Q -6 1Q -5 

Flux (>100MeV) [photons cm" 2 sec" 1 ] 

Fig. 1. — The upper panel: the cumulative flux distribution of gamma-ray emitting clusters of galaxies. The three curves are the theoretical 
predictions using three different formulae for the formation rate of dark haloes (-Rform)i i-e., the derivative of the Press-Schechter (PS) mass 
function (solid line), formation rate of Sasaki (1994; dot-dashed line) and Kitayama & Suto (1996a; dashed line). The dotted line is the 
flux distribution of uniformly distributed sources in the Euclidean space. The observed distribution of the unidentified EGRET sources with 
\b\ > 45° is shown by the thick (only the 'em' sources) and thin (all sources including non-'em') lines. The lower panel: the solid, dashed, and 
dot-dashed lines are mean redshift, mass [10 16 Mq], and angular radius [degree] of gamma-ray clusters brighter than a given flux, respectively. 
The dotted line is the contribution of gamma-ray clusters brighter than a given flux to the diffuse extragalactic gamma-ray background at 
100 MeV. The derivative of the PS mass function is used. The sensitivity limits of the EGRET and GLAST experiments are shown in the 
figure. 



ber of 'steady' unidentified sources with |6| > 45° defined 
by Gehrels et al. (2000). This result is robust against 
changes in the adopted prescription for i?f or m- In the lower 
panel of Fig. [t], we also show mean mass, redshift, and an- 
gular radius # v ir corresponding to r v ; r of such gamma-ray 
clusters brighter than a given flux. These quantities are 
M ~ 1O 15 M , z ~ 0.05, and 8 vil ~ 1° for clusters above 
the EGRET sensitivity limit. Considering the EGRET 
angular resolution, the typical radius of ~ 1° is consis- 
tent with the fact that a significant fraction of isotropic 
unidentified sources are indicated as possibly extended. 

It is predicted that more than a few thousands of form- 
ing clusters will be detected by future missions such as 
the GLAST (Gehrels & Michelson 1999), and the im- 
proved angular resolution may reveal the extended profile 
for nearby gamma-ray clusters with higher statistical sig- 
nificance. Another important prediction is that GLAST 
will observe the flattening of the log AT-log F curve due to 
the cosmological effects, compared with the expectation 
of a uniform source distribution in the Euclidean space 
(dotted line in the upper panel of Fig. [|). 



4. EXTRAGALACTIC GAMMA-RAY BACKGROUND 

Our formulation also allows us to calculate the EGRB 
flux and spectrum as 

—l= dz dM 7V ?' ' R lorm (M,z) — (10) 
de 7 J J ae 7 dz 

where (dn 7 /de 7 ) is the gamma-ray number density that is 
related to the EGRB flux as (dF/dej) = c((in 7 /c?e 7 )/(47r), 
and t is the cosmic time. This flux becomes £j(dF '/ 'de 7 ) <~ 
f keV cm~ 2 sec" 1 sr _1 at f 00 MeV, in good agreement with 
the observation (Sreekumar et al. f 998) as well as the pre- 
vious simpler estimation assuming that the average tem- 
perature of baryons in the universe is ~ keV at present 
(Loeb & Waxman 2000). In fact, we have checked that 
the mass-averaged temperature of virialized haloes in the 
universe as a function redshift, that is calculated by the 
PS theory, agrees within a factor of 2 with a numerical 
simulation (Cen & Ostriker 1999) on which the previous 
EGRB estimate was based. This fact gives a justification 
for the use of the PS theory to calculate the gamma-ray 
emitting objects. The dotted line in the lower panel of 
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Fig. 2. — The spectrum of the diffuse extragalactic gamma-ray background radiation. The data is from Sreekumar et al. (1998). The 
parameter £g is the fraction of magnetic energy in the gravitational energy given to baryonic gas in a collapsed halo, with £b = 10 _1 
(dashed line), 10 — 3 (solid), and 10 -5 (dot-dashed). The absorption of gamma-rays above ~ 100 GeV by the pair-production interaction with 
the cosmic infrared background is not taken into account in the upper panel, while it is in the lower panel. (The secondary gamma-rays 
reprocessed by the produced pairs are not taken into account in either of the panels, see text.) 



Fig. shows the contribution to the EGRB at 100 MeV 
by objects brighter than a given flux. We predict that 
GLAST will resolve about 20-30 % of the EGRB as dis- 
crete gamma-ray clusters, if structure formation is the ma- 
jor origin of the EGRB. 

The strength of magnetic field in the shocked baryons 
is important for the question whether the EGRB spec- 
trum extends up to ~ 100 GeV as observed. In Fig. |^ 
we show the EGRB spectrum with several values of £g. 
This result shows that the magnetic field strength corre- 
sponding to £b ~ 10~ 5 of the baryon energy density is 
sufficient for the EGRB to extend beyond 100 GeV. The 
magnetic field observed in the intracluster gas (£b ~ 10~ 3 ) 
is much stronger than this, and it is also theoretically rea- 
sonable to expect that the turbulent motion in collapsed 
objects amplifies the seed magnetic field made by the bat- 
tery mechanism well beyond £s > 10~ 5 within the dy- 
namical time (Kulsrud et al. 1997). Therefore, physically 
reasonable magnetic field strength can explain the exten- 
sion of the EGRB spectrum beyond 100 GeV, and it is 
likely that a considerable fraction of gamma-rays above 
100 GeV is absorbed by the interaction with the cosmic 
infrared background radiation producing electron-positron 
pairs. The effect of intergalactic absorption is shown in 



the lower panel of Fig g, using the optical depth of inter- 
galactic pair-production presented in Totani (2000). The 
absorbed TeV gamma-rays will be reprocessed into GeV 
gamma-rays by the produced pairs, and distort the EGRB 
spectrum (e.g., Coppi & Aharonian 1997). Although these 
secondary GeV gamma-rays are not taken into account 
here and it is beyond the scope of this paper, it is impor- 
tant to study how large is this spectral distortion in future 
work. 

5. DISCUSSION 

We here discuss the expected properties of gamma-ray 
clusters of galaxies. Perhaps the most natural question in 
this regard would be "Are they already observed in other 
wavebands such as x-rays or optical surveys?" We have 
checked that there is no statistically significant association 
of the ROSAT Brightest Cluster Sample (RBCS, Ebeling 
et al. 1998), within the 95% error circles of the unidentified 
sources with \b\ > 30° in the EGRET catalog. We have 
also checked the correlation with the clusters in the re- 
vised Abell catalog (Abell, Corwin, & Olowin 1989), and 
no statistically significant associations are found, either. 
However, in the following we argue that the gamma-ray 
clusters proposed in this paper are very difficult to detect 
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in x-rays or optical bands compared with ordinary clus- 
ters identified in these wavebands, and hence our scenario 
is not rejected by these results. 

5.1. Detectability of gamma- ray clusters in X-rays 

We first estimate the expected x-ray flux from gamma- 
ray clusters. Baryonic gas in most clusters of galaxies 
observed in x-rays seems to be in approximate hydro- 
static equilibrium with the surface brightness well fit- 
ted by a density profile, p gas (r) oc [1 + (r/r c ) 2 ] _1 (e.g., 
Sarazin 1988), where r c is the core radius that is typ- 
ically about ~ 10 times smaller than the virial radius. 
Since the x-ray emissivity is proportional to /o gas , the x- 
ray emission is strongly concentrated into the central re- 
gion. Assuming the above density profile and the self- 
similar model as described in Kitayama & Suto (1996b), 
1 a typical gamma-ray cluster detectable by the EGRET 
with M ~ 10 15 M Q and z ~ 0.05 would have the x-ray flux 
2.4 x lO" 11 erg cm~ 2 s" 1 in 0.1-2.4 keV. 

The inverse-Compton flux is also expected to be com- 
parable with the thermal emission. By equating tjc and 
ishock in §H we get the cooling photon energy e 7iCOO i = 
2.0(1 + z)~ 5 keV, below which the electron cooling time 
is longer than the dynamical time. Then the IC spectrum 
extends down to around x-ray band with dN 1 /dt 1 oc e^ 2 , 
while it becomes harder at wavelengths longer than x-rays 
with d/V 7 /de 7 oc e~ 15 . If the gamma-ray flux at 100 MeV 
is ~ 10 -7 photons cm" 2 s _1 that is the EGRET threshold, 
the IC x-ray flux {vF v ) is ~ 1.6 x 10 -11 erg cm -2 s _1 . 
Therefore, the thermal and IC fluxes are well above the 
flux limit - 4 x 10^ 12 erg cm~ 2 s" 1 of the RBCS. 

However, it takes nearly the dynamical time for the clus- 
ter gas to reach hydrostatic equilibrium after the collapse, 
and gamma-rays from the shock generated by the gravita- 
tional collapse are radiated away within that period. Then 
it is likely that the density profile of gamma-ray emitting 
clusters is more irregular and extended than ordinary x- 
ray clusters. In fact, if the unidentified 'em' sources in the 
EGRET catalog are actually extended, they must have 
typical angular size of about degree, from the source loca- 
tion accuracy of the EGRET. As we have shown, angular 
size of about 1° is theoretically reasonable if the emission 
is extended to the virial radius. When the density pro- 
file is not concentrated into the central region but rather 
constant within the virial radius, the x-ray luminosity be- 
comes lower than the self-similar model by a factor of ~ 3.7 
because of the lower central density. 

Furthermore, the surface brightness of such loose clus- 
ters should be drastically dimmer than ordinary x-ray clus- 
ters. In the self-similar model with r v ; r 3> r c , the core 
gas density is p gaS:C ~ (l/3)(r vir /r c ) 2 p gaS:Vir - 50p gaS:Vir , 
where p gas ,vir = (^s/^o)Pvir is the virial gas density that 
is the average gas density within r v ; r . On the other hand, 
if the gas density profile of gamma-ray clusters is roughly 

1 Here we have assumed that the core radius is proportional to the 
with r v ; r (10 15 MQ, 0) = 2.6Mpc, where the normalization is chosen to m 
1984). 



constant at p gas , v ir out to r v ; r , the x-ray surface bright- 
ness of such a loose cluster is dimmer than the central 
surface brightness of the self-similar model by a factor of 

~ (r c Avir)(/°gas,c/Pgas,vir) 2 ~ 2 °0: s i nce tne x-ray emis- 
sivity is proportional to p gas . It crucially affects the de- 
tectability of x-rays from gamma-ray clusters. The de- 
tectability of x-rays should be described by the signal- 
to-noise ratio (S/N) against the x-ray background flux. 
The noise level is proportional to (image area) 1 / 2 , and 
hence S/N oc F/r, where F and r are the flux and the 
image radius, respectively. We have compared the value 
of F/r of the extended gamma-ray clusters detectable by 
the EGRET and those of the clusters in the RBCS. We 
found that the F/r of gamma-ray clusters is by a factor 
of 3 smaller than the minimum F/r of the RBCS clusters. 
The absence of association between the RBCS and the 
EGRET sources is therefore not in contradiction to our 
scenario. On the other hand, deeper observation of candi- 
date gamma-ray clusters by Newton, for example, might 
detect the x-ray emission extended to about 1° with the 
flux estimated above, that would provide a clear test of our 
scenario. Such x-ray emission should reflect the structure 
of shocks in dynamically forming clusters, and imaging 
study is of great interest. 

5.2. Detectability in the Optical Surveys 

Here we again emphasize that the gamma-ray clusters 
are expected to be more extended than clusters that have 
already stabilized. It is known that the surface density 
profile of galaxies in a cluster can well be described by the 
King profile, a(r) oc [1 + (r/rc) 2 ]^ 1 with the core radius of 
^100 kpc that is comparable with the core radius of x-ray 
profile (e.g., Adami et al. 1998). If we assume a roughly 
constant surface density out to ~ r v j r rather than the King 
profile for gamma-ray clusters, the average surface density 
Cav ~ -^gai/(""' , vir) i s lower than the central surface density 
of the King profile a c ~ N ga i/[2nr 2 ln(r v ; r /r c )], by a factor 
of cr c /cr av ~ [2\n(r vil ./r c )}~ 1 (r vil /r c ) 2 - 30. Here iV ga i is 
the total number of galaxies within r v ; r . This dimming 
factor is not so significant as that for the x-ray surface 
brightness, but that should make the optical identifica- 
tion very difficult because of the contamination by fore- 
ground and/or background field galaxies. Therefore, we 
consider that no statistically significant association with 
the known optically identified clusters does not immedi- 
ately reject our scenario. Instead, it is necessary to study 
in the future the correlation between the EGRET sources 
and galaxy catalogs taking into account the possibility 
that the gamma-ray clusters are considerably extended. 
Search in optical bands has an advantage over the search 
in x-rays, in a sense that the dimming of surface num- 
ber density compared with ordinary clusters is less severe 
than x-rays whose emissivity is proportional to p gas . The 
typical density of such loose clusters is close to the virial 

virial radius as r c = 0.21(h/0.7) _1 Mpc X [r vlr (M, z)/r vil (10 15 Mq, 0)] 
atch the local observations (Abramopoulos & Ku 1983; Jones & Forman 
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density, that is about a few hundreds times higher than 
the mean density of the universe. 

We have also noticed that there are a considerable num- 
ber of 'em' sources in the EGRET sources identified as 
AGNs. If they were actually extended sources, it might be 
speculated that some of them are also gamma-ray clusters 
including an AGN as a member galaxy. Time variability 
of these sources would be an important test to check this 
possibility. 

5.3. On the recent follow-up observations for 3EG 
J1835+5918 

Recent follow-up observations by Mirabal et al. (2000) 
for one of the high-latitude unidentified EGRET sources 
(3EG J1835+5918) have found a diffuse x-ray emission 
from an uncatalogued cluster of galaxies at z = 0.102. Al- 
though this cluster is outside the 99% error ellipse of 3EG 
J1835+5918 whose radius is 12', the separation between 
the centers of the x-ray cluster and 3EG J1835+5918 is 
about 0.65°, that is within our expectation of the typical 
angular radius of gamma-ray clusters detectable by the 
EGRET, ~ 1°. As discussed above, x-ray emission is ex- 
pected from a region where the intracluster gas reached 
hydrodynamical equillibrium, while gamma-rays are emit- 
ted from a region still hydrodynamically unstable. There- 
fore, it is not surprising that the positions of x-ray and 
gamma-ray emissions are different unless the separation is 
well beyond the virial radius of ~ 1°. 

3EG J1835+5918 is not an 'em' source, and this source 
may not be an extended source. The radius of the 99% 
confidence ellipse, 0.2°, is considerably smaller than the 
expected angular size of gamma-ray clusters. However, 
high energy electrons emitting GeV gamma-rays have very 
short life time compared with the shock propagation time 
(tic ~ 10~ 3 £ s hock), and gamma-ray emitting region may 
be very clumpy in a cluster. (On the other hand, x-ray 
emitting electrons have a cooling time comparable with 
or longer than i s hock for IC and thermal radiations, and 
hence x-ray emitting region should be much less clumpy 
and extended with the size ~ 9 v i T , as discussed in § |5.1| .) 
Therefore, it is possible that the gamma-ray size of the 
3EG J1835+5918 is considerably smaller than the phys- 
ical size of a whole forming cluster. This consideration 
also suggests a possibility that some of gamma-ray clus- 
ters may be observed as multiple sources within ~ 9 v i T , 
that may be revealed by future gamma-ray missions. 

One of the characteristics that the source of 3EG 
J1835+5918 must have is very weak radio flux that is at 
least two orders of magnitudes fainter than any of the se- 
curely identified EGRET blazars (Mirabal et al. 2000). 
The spectrum of blazars is well understood by the two 
components of radiation by the same population of non- 
thermal electrons, i.e., synchrotron radiation in radio, op- 
tical, and X-ray bands, and inverse-Compton radiation 
in GeV and TeV gamma-ray bands (e.g., Inoue & Taka- 
hara 1996; Kataoka et al. 1999). The ratio of luminosi- 



ties by the two processes is, as is well-known, given by 
the ratio of magnetic energy density to the target pho- 
ton energy density, that is typically of order unity for 
blazars. On the other hand, this ratio is Ub/Ucmb ~ 
2.7 x 10- 3 (£ s /10- 3 )(M/10 15 Af Q ) 2 / 3 for gamma-ray clus- 
ters, by using the expression of B given in §||. Therefore, 
the Ub/Ucmb ratio is generally much smaller than the 
unity, and very weak radio flux compared with identified 
blazars can be reasonably explained. 

Based on the above arguments, we suggest that the un- 
catalogued x-ray cluster near 3EG J1835+5918 may be a 
gamma-ray cluster proposed in this paper, and further ob- 
servations for this cluster and the surrounding region are 
very important. Our model predicts that a cluster emit- 
ting a flux of ~ 6.06 x 10~ 7 photons cm _2 s _1 above 100 
MeV (the flux of 3EG J1835+5918, Hartman et al. 1999) 
at z = 0.102 should have a total mass of ~ 7 x 1O 15 M 
and r vir ~ 5 Mpc (6> vir - 0.7°). 

6. CONCLUSIONS 

In this paper we have proposed a new candidate of 
unidentified EGRET sources: gamma-ray clusters that 
are just dynamically forming and emit gamma-rays due 
to inverse-Compton scattering of CMB photons by shock- 
accelerated electrons. Based on the standard theory of 
structure formation and assuming the injection of ~5% of 
shock energy at the formation into nonthermal electrons, 
we have shown that a few tens of such clusters should have 
already been detected by EGERT, and a significant frac- 
tion of the isotropic component of unidentified EGRET 
sources can be accounted for. Such gamma-ray clusters 
are expected to be very extended; the x-ray surface bright- 
ness and surface number density of galaxies could be lower 
than those of ordinary clusters by a factor of ~ 200 and 
~ 30, respectively. Therefore it should have been very dif- 
ficult to detect gamma-ray clusters in the past x-ray or 
optical surveys, and our scenario is in accord with appar- 
ent no-associations between unidentified EGRET sources 
and x-ray or optical clusters. 

It will be of great significance to perform x-ray or optical 
observations to search for such loose clusters of galaxies in 
the regions of high-latitude unidentified EGRET sources. 
The future gamma-ray projects such as GLAST will also 
provide a direct test of our scenario. If our scenario is 
true, a new population of "gamma-ray clusters" will pro- 
vide us in the future a new probe of dynamically evolving 
structures in the universe that cannot be traced by x-ray 
or optical clusters of galaxies. 
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